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Abstract With the growing efforts in isolating solid-state qubits from external decoherence sources, the material-inherent 
sources of noise start to play crucial role. One representative example is electron traps in the device material 
or substrate. Electrons can tunnel or hop between a charged and an empty trap, or between a trap and a gate 
electrode. A single trap typically produces telegraph noise and can hence be modeled as a bistable fluctuator. 
Since the distribution of hopping rates is exponentially broad, many traps produce flicker-noise with spectrum 
close to 1//. Here we develop a theory of decoherence of a qubit in the environment consisting of two-state 
fluctuators, which experience transitions between their states induced by interaction with thermal bath. Due 
to interaction with the qubit the fluctuators produce 1 //-noise in the qubit's eigenfrequency. We calculate the 
results of qubit manipulations - free induction and echo signals - in such environment. The main problem is that 
in many important cases the relevant random process is both non-Markovian and non-Gaussian. Consequently 
the results in general cannot be represented by pair correlation function of the qubit eigenfrequency fluctuations. 
Our calculations are based on analysis of the density matrix of the qubit using methods developed for stochastic 
differential equations. The proper generating functional is then averaged over different fluctuators using the 
so-called Holtsmark procedure. The analytical results are compared with simulations allowing checking accu- 
racy of the averaging procedure and evaluating mesoscopic fluctuations. The results allow understanding some 
observed features of the echo decay in Josephson qubits. 

Keywords: Qubits, Decoherence, 1 //-noise 



1. Introduction and model 



The dynamics of quantum two-level systems has recently attracted special attention in connection 
with ideas of quantum computation. A crucial requirement is to the phase coherence in the presence 
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of noisy environment [1]. Solid state devices have many advantages for realization of quantum compu- 
tation that has been confirmed by several successful experiments, for a review see, e. g., Ref. [2] and 
references therein. In solid-state realizations of quantum bits (qubits) the major intrinsic noise is due to 
material-specific fluctuations (substrate, etc). Concrete mechanisms of these fluctuations depend upon 
the realization. In particular, in the case of charge qubits the background charge fluctuations with 1// 
spectrum are considered as most important [2]. They are usually attributed to random motion of charges 
either between localized impurity states, or between localized impurity states and metallic electrodes. 

The conventional way to allow for the noisy environment is to describe it as a set of harmonics oscil- 
lators with a certain frequency spectrum. The resulting "spin-boson models" were extensively discussed 
in the literature, see for a review Refs. [3] and [4]. Applications of these models to concrete qubit 
implementations have been recently reviewed by Shnirman et al. in Ref. [5]. 

1.1 Spin-boson model 

Conventionally, the quantum system which we will call the qubit is assumed to be coupled linearly to 
an oscillator bath with interaction Hamiltonian 



H = a z X , X = Y,Cj (bj + St) . (1.1) 



Here o~i, i = x,y, z, are the Pauli matrices describing the qubit, while bj and Sj stand for bosons. The 
decoherence is then expressed in terms of the symmetric correlation function 

S x (u) = (jx(t),X(0)] + J = 2 J(u) coth ^ , (1.2) 

where J(uj) the bath spectral density, J(uS) = tt J2j Cf 5{lo — ujj). Here and below we put h = 1 and 
/cb = 1- In the simplest case the decoherence can be characterized by [5] 



K{t) = - InTr r e -<*(«) e <*(o) pA (1.3) 

where pb is the density matrix of the thermal bath, and the bath phase operator is defined as 

Ht) = ^(2CjM) e iHbt (St - bj) e~ iHbt , (1.4) 
3 

Jib being the bath Hamiltonian. The quantity K,(t) is conventionally expressed through the bath spectral 
density, J(uj), as [3] 
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The most popular assumption about J(u), namely J(oj) = (o>ttuj/2)@(lj c — uS), is "Ohmic dissipation". 
Here a is a dimensionless coupling strength, while <d(x) is the Heaviside unit step function. 

Important features of this approach is that (i) the decoherence is determined solely by the pair corre- 
lation function Sx(lo) that assumes the noise to be Gaussian; (ii) Sx(u>) is related to the bath spectral 
density through the fluctuation-dissipation theorem, which assumes the system to be equilibrium. As 
long as these assumptions hold, the method provides powerful tools to analyze the decoherence. 

Several attempts, see Ref. [5] and references therein, were made to extend the spin-boson model to 
the so-called "sub-Ohmic" case, in particular, to the case of l//-noise where Sx(w) °c \oj\~ . 

We believe that 1/ f noise is a typical nonequilibrium phenomenon. It is due to the fact that some 
excitations of the environment relax so slowly that cannot reach the equilibrium during the measuring 
time. As a result, the fluctuation-dissipation theorem cannot be applied. Moreover, 1// noise is not a 
stationary Markov process. Indeed, it is created by the fluctuators with exponential broad distribution of 
the relaxation rates and thus is not fully characterized by its pair correlation function. 
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Figure 1.1. Schematic sketch of the distribution of charged traps. They are located near the gate's surface and produce 
oppositely charged images. 



1.2 Spin-fluctuator model 

Several attempts were made to study the role of non-Gaussian and non-Markovian nature of the 1//- 
noise for various examples of coherent quantum transport such as resonant tunneling [6], ballistic trans- 
port through a quantum point contact [7], Josephson effect [8], Andreev interferometer [9]. In connec- 
tion to qubits, a similar model has been recently studied by Paladino et al. [10]. In this paper dynamical 
charged traps were considered as two-level systems (TLS) with exponentially-broad distribution of hop- 
ping rates, and the "free induction signal" of a qubit was numerically analyzed for a narrow distribution of 
the coupling constants, Vi, between the traps and the qubit. Quantum aspects of non-Markovian kinetics 
were addressed in Ref. [11]. 

The aim of the present work is to revisit this problem. In the following we will consider a spin- 
fluctuator model, similar to that considered in Ref. [10], which takes into account both nonequilibrium 
and non-Markovian effects. Analysis of this model shows that nonequilibrium effects are important. In 
particular, we will address the role of the distribution of coupling constants between the fluctuators 
and the qubit. The distribution of is probably broad for most possible devices and realistic situations 
because the fluctuators are located at different distances from the qubit. We will show that a broad 
distribution of the coupling constants leads to significant modifications of the decoherence dynamics. 

The broad distribution of the coupling constants, makes the model, which we will consider, essentially 
similar to the conventional models of the spectral diffusion in glasses. The concept of spectral diffusion 
was introduced by Klauder and Anderson [12] as early as in 1962 for the problem of spin resonance. 
They considered spins resonant to the external microwave field (spins A) which generate echo signals, 
and surrounding non-resonant spins (spins B). Due to interaction between spins A and B, stochastic 
flip-flops of spins B lead to a random walk of A-spins frequencies. This random walk is referred to as 
the spectral diffusion. 

Black and Halperin [13] applied the concept of spectral diffusion to low-temperature physics of 
glasses. They used ideas of Ref. [12] to consider phonon echo and saturation of sound attenuation 
by two-level systems [14] in glasses. Important generalizations of these results were made by Hu and 
Walker [15], Mainard et al. [16] and Laikhtman [17]. 

Following these ideas, we assume that the qubit is a two-level system (TLS) surrounded by so-called 
fluctuators, which also are systems with two (meta)stable states. One can imagine several realizations 
of two-level fluctuators. Consider, e. g., structural dynamic defects, which usually accompany really 
quenched disorder. These defects, if not charged, behave as elastic dipoles, i. e., they interact with 
the qubit via deformational potential. The interaction strength in this case depends on the distance, r, 
between the fluctuator and the qubit as r~ 3 . Consequently, the distribution of the coupling constants is 
V(v) oc t>~ 2 . Charge traps near the gates also produce dipole electric fields, see Fig. 1.1. In this case 
v = e 2 (a ■ r)/r 3 oc r~ 2 . Assuming that a <sCr (otherwise charge would tunnel directly to the qubit and 
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the device will not work) and integrating the fluctuators contributions along the 2D gate surface we again 
obtain V(v ) oc v~ 2 . It is exactly the distribution of the coupling constants in glasses, where two-level 
systems (TLS) interact via dipole-dipole interaction [13]. 

It is crucial that due to their interaction with the environment the fluctuators switch between their 
states. This switching makes the fields acting upon the qubit time-dependent. The decoherence is caused 
by the time-dependence of the random field - a static field can only renormalize the qubit's interlevel 
spacing. The dynamics of the fluctuators is relaxational, so they are rather "relaxators" than oscillators. 
The decoherence is thus determined by the fluctuators relaxation rates, which should be compared to the 
measurement time. At low temperatures the fluctuators are frozen in the ground states, and their dynam- 
ics is slow. Therefore, the fluctuators-induced decoherence significantly decreases with the temperature. 

The paper is organized as follows. In Sec. 1.1.3 we reformulate the model for spectral diffusion 
in glasses for the case of a qubit and consider a simplified version of the theory. This approach is, 
strictly speaking, applicable only to the fluctuators with the small interlevel spacings, E <^ T, and 
the numerical factors that follow from this approximation should not be trusted. However the resulting 
qualitative physical picture is believed to be correct, because the fluctuators with E > T are frozen in 
their ground states, they thus do not fluctuate and do not contribute to the decoherence. The problem of 
a single neighboring fluctuator is considered in Sec. 1.2 followed by the discussion of the averaging over 
the ensemble of fluctuators in Sec. 1.3. In Sec. 1.5 we show that different fluctuators are responsible for 
the decoherence and for the flicker noise. Consequently, the decoherence cannot in general be expressed 
through the pair correlation function of random fields acting upon the qubit. 

1.3 Detailed description 

A qubit coupled to the environment will be modeled by the Hamiltonian 

H = H q + H mw + H q F + TCf , (1.6) 

where H q and Hf describe the qubit and the fluctuators separately. A completely isolated qubit is just 
a system that can be in one of two states and is characterized by the energies of these states and the 
tunneling probability between them. ~H q can thus be written as the Hamiltonian of the qubit pseudospin 
in a static "magnetic field" B = {B x , B z }, where B z characterizes the splitting of the energies of the 
two states, and B x is responsible for the tunneling. (Parallel shift of the qubit energies is, of course, 
irrelevant). One can diagonalize such a Hamiltonian, by simply choosing the direction of the new z-axis 
to be parallel to B and write the rotated Hamiltonian TL q as 

H q = ~Ba z . (1.7) 

Tip in turn can be diagonalized and split into three parts 

Jip = Hp + Hp-em + H env . (1.8) 
We model fluctuators by two-level tunneling systems, which Hamiltonians can also be diagonalized: 

nf = \Y,Eiri (1.9) 

i 

where the Pauli matrices tW correspond to z-th fluctuator. The spacing between the two levels, E{ is 
formed by the diagonal splitting, Aj, and the tunneling overlap integral, Aj 

Ei = ^Af + A? = Ai/sinfli. (1.10) 

To account for the flip-flops of the fluctuators one needs to include environment. To be specific we 
model the environment by a bosonic bath. This applies not only to phonons, but also to electron-hole pairs 
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in conducting part of the system [18]. Neglecting the interactions of the fluctuators with the environment 
that do not cause the flip-flops we specify the environment-related parts of the Hamiltonian as 

n em = 5>m(S)& + 1/2) , (lid 

i n 

where \i represents the boson quantum numbers (wave vector, etc.), and ^ are the constants of the 
coupling between the fluctuators and the bosons. 

It is crucial to specify the interaction, 7i q F, between the qubit and the fluctuators. Following Ref. [13], 
we assume that 

H q F = ^ Vj cr z rl , Vi = g(ri)A(iii) cos 8i. (1.13) 

i 

Here 9i is determined by Eq. (1.10), rij is the direction of elastic or electric dipole moment of ith fluc- 
tuator, and n is the distance between the qubit and ith fluctuator. The functions A(rij) and g{ri) are not 
universal. 

Interaction between the fluctuators and the environment manifests itself through time-dependent ran- 
dom fields £i(t) applied to the qubit. These low frequency, uj <C T, fields can thus be treated classi- 
cally. Accordingly, one can substitute Tt q F by the interaction of the qubit pseudospin with a random, 
time-dependent magnetic field, X\{t), which is formed by independent contributions of surrounding 
fluctuators: 

U qF = Xi(t) a z , Ai(t) = ^viti(t) , (1.14) 

i 

so that Vi, Eq. (1.13) determines the coupling strength of the i-th fluctuator with the qubit, while the 
random function characterizes the state of this fluctuator at each moment of time. Below we assume 
that fluctuator switching itself is an abrupt process that takes negligible time, and thus at any given t either 
&(t) =0ar&(t) = 1. 

Note that we directed X\(t) along z-axis, i.e., neglected possible transitions between the qubit eigen- 
states induced by the random fields. This can be justified by the low frequency of these field. (From 
the practical point of view a qubit that switches frequently by the external noise does not make a decent 
device.) 

To manipulate the qubit one should be able to apply ac "magnetic field", F(t). In general this field 
not only causes the transitions between the qubit eigenstates (F x ), but also modulates the qubit level 
spacing B in time (F z ). The latter effect is parasitic and should be reduced. Here we simply neglect this 
modulation and assume that F is parallel to the x-axis (F x = F). Accordingly 

Hr Dm = {l/2)F(t)-a x . (1.15) 

For the manipulation to be resonant, the frequency Q of the external field F should be close to B. 
Combining Eqs. (1.7), (1.9), (1.14), and (1.15) we substitute the initial Hamiltonian (1.6) by: 

H = \ [Eq + X(t)] a z + X -F(t)a x + ^ E iT \ . (1.16) 

i 

Here Eq determines the original (t = 0) value of the eigenfrequency of the qubit, while X(t) is the 
random modulation of this eigenfrequency caused by the flips of the fluctuators: 

Eq = B + X 1 (0), X(t) = X±(t) — Xi(0) . (1.17) 

We have not explicitly included into the Hamiltonian (1.16) the interaction (1.11) of the fluctuators 
with phonons or electrons, which causes flips of the fluctuator. We account for these flips by introducing 
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finite transition rates between the fluctuator's states. The transition rates can be calculated in the second 
order of the perturbation theory for the fluctuator-phonon/electron interaction [19, 18]. As a result, the 
ratio, 7 + /7~, of the rates for the upper, 7 + , and the lower, 7", states in the equilibrium is determined 
by the energy spacing between these two states: 7^/7" = exp (Ei/T). Accordingly, these rates can be 
parameterized as 

7± = (l/2) 7i (T, Ei) [1 ± tanh^/T)] sin 2 { (1.18) 

where 0j is given by Eq. (1.10). 

The dependence of 7o(-E'i) on Ei is determined by the concrete relaxation mechanism: for phonons [19] 
7o(-Ei) oc Ef, while for electrons [18] jo(Ei) oc E{. Note that the average value of £j(t) depends only 
on Ei and T: 

(£(*)>*-*» = 7~/(7 + + 7") = [1 + exp (Ei/T)]- 1 . (1.19) 

There can also be a direct interaction of the qubit with the bosonic bath. One can introduce the transition 
rate 7 9 (T) due to this interaction in the same way as 7 i(T). 

Below we will often use a simplified version of the theory assuming that the only fluctuators that 
contribute to dephasing are those with Ei <C T. In this approximation 

jf = (1/2)70 (T) sin 2 6i, (1.20) 

70 is thus the maximal fluctuator relaxation rate at a given temperature T. This assumption significantly 
simplifies the formulas and produces, as one can show [17, 20], correct order-of-magnitude results. 

The model formulated above differs from the spin-boson models by statistics of the random force 
X(t). It allows one to describe the qubit decoherence in a simple way without loosing track of the 
essential physical picture. 

1.4 Qubit manipulations 

We parameterize the qubit's density matrix as 



n —if e 



iStt 



if* e~ int 1 - n 



(1.21) 



and commute it with the Hamiltonian, Eq. (1.16), using the resonance approximation. (The frequency, Q 
of the applied field F is assumed to be close to the qubit eigenfrequency, Eq.) Including also the inherent 
qubit relaxation (7^) we obtain the following equations of motion: 

= -2 7q (n-n )-FRef, (1.22) 



Of 
dj_ 

dt 



i [Eq + X(t) -n}f- lq f + F (n-^j. (1.23) 



These equations have been obtained in Refs. [21, 22] for the problem of spectral diffusion in glasses, F 
in Eq. (1.16) playing the role of the Rabi frequency of the resonant pair. 

The equation set (1.22) belongs to the class of stochastic differential equations. In the following we 
use the methods [23-29] developed for these equations to study the qubit response to various types of 
the manipulation. 

Currently the experimentally observable signal is an accumulated result of numerous repetitions of the 
same sequence of inputs (e. g., pulses of the external field). To be compared with such measurements, 
solutions of Eq. (1.21) should be averaged over realizations of the stochastic dynamics of the fluctuators. 
Provided that the time intervals between the sequences of inputs are much longer than the single-shot 
measuring time we can separate the averaging over the initial states of the fluctuators, £j(0), from the 
averaging over their stochastic dynamics, X(t). 
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In the absence of the external ac field n(t) should approach its equilibrium value given by the Fermi 
function, uq(Eq) = [1 + exp(Eo/T)]~ 1 , while the off-diagonal matrix element of the density matrix 
should tend to zero. If the external ac field is switched off at t = then the solution of Eq. (1.21) has the 
form 

n(t) = n (E ) + [n(0) - noiEo)^- 2 ^ , (1.24) 
f(t) = f{{S) e ~^ t+i{ ~ E °- n)t+i ^ x(t')dt' _ ( 12 5) 

We need to average this signal, known as the free induction signal over both Eq and X(t). If the time 
t after the pulse is short enough, then the fluctuators remain in their original states and X(t) can be 
neglected. At t = the system of the fluctuators is supposed to be in the equilibrium, and thus probability 
for £i = 1 is no(Ei). Substituting Eq from Eq. (1.17) we find that exp[i(Eo — Q)t] averaged over the 
initial realization equals to 

^ e i(E -Q)t\ = e i(B-a)t JJ [i _ no (Ej) + n Q {E j )e iv i t ] . (1.26) 

j 

It follows from Eq. (1.26) that the observed free induction signal involves the oscillations with fre- 
quencies that differ from B — by various combinations of Vj. As a result, in the presence of a large 
number of the fluctuators the free inductance signal decays in time even when X(t) = 0, i. e. when 
the fluctuators do not switch during one experimental run. However, this decay has little to do with the 
decoherence due to irreversible processes. 

Much more informative for studies of genuine decoherence are echo experiments, when the system 
is subject to two (or three) short ac pulses with different durations t\ and 72, the time interval between 
them being t\2 (or t\i and 713, respectively). Considering echo, we assume that the external pulses are 
short enough for both relaxation and spectral diffusion during each of the pulses to be neglected. The 
echo decay is known to be proportional to the "phase-memory functional" [30] 

¥0(f), t) = ^exp (i J* 0(t')X(t') dt'^j ^ . (1.27) 

Here for 2-pulse echo t = 2t\2 and 

{0 for t' < , 
1 for < t' < ri2 , (1-28) 
-1 for r 12 <t'. 

In the case of the 3-pulse echo one would put t = t\% + T13 and 

m = k 



The functional (1.27) can be used to describe the decoherence of the free induction signal, substitut- 
ing 0(f) = 0(0- in this case, however, one should understand X(t') as Yl,i v i^i^') rat her than use 
Eq. (1.17). For the echo experiments t in Eq. (1.27) is chosen in such a way that the integral of 0{t') 
from zero to t vanishes. As a result, the time-independent part of X, i. e., dispersion of £j(0) becomes 
irrelevant. Below we evaluate the phase-memory functional (1.27) for the free induction as well as for 
schemes of the measurement. 






for 


t' < 0, 


1 


for 


< t' < T12 , 





for 


7l 2 <t' < 713 


-1 


for 


T13 < t' . 



2. Results for a single fluctuator 

2.1 Random telegraph noise 

The process, which is described by a random function, £(t), that acquires only two values: either 
£ = or £ = 1 is known as random telegraph process. In this section we assume that Ei 2> T and thus 
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the two states of each fluctuate* are statistically equivalent, i.e., the time-average value of £(£) equals to 
<£(i))^oo = 1/2. 

To evaluate the memory functional (1.27) we first introduce auxiliary random telegraph processes 
defined as 

Z± (t) = ±(-ir^\ (i.30) 

where n(ti, t 2 ) is a random sequence of integers describing number of 'flips" during the period {t\,t 2 ), 
so that n(t,t) = 0. The fact that = 1 substantially simplifies the calculations. The 'flips" of a 

given fluctuator induced by its interaction with the bosonic bath should not be correlated with each other. 
Accordingly, n(t\, t 2 ) obeys the Poisson distribution, i. e., the probability, V n (ti,t2)), that n(t\, t 2 ) = n 
equals to 

VnQtu h)) = (n(tl ', t2))ra e-<"frA» , (n(t h t 2 )) = 7 \ tl - t 2 \ , 



nl 



(1.31) 



where 7 has a meaning of the average frequency of "flips". From Eqs. (1.31) and (1.30) it follows that 

(z±(t)) = = ±e- 2 ^l , 



n=0 



e -2 7 ( tl - ta ) j h > t2 



<*fc(*l)*fc(k)> = ((-lf^) 

It is convenient to describe different measurement schemes by making use of the generating function- 



al 



i>±\P,t] = (^exp 



(1.32) 



z ± (t) 



where f3(t') is the same function as in Eqs. (1.27), (1.28), and (1.29), while the constant v will later play 
the role of the qubit-fluctuator coupling constant. To evaluate the functionals (1.32) consider a set of the 
correlation functions 

M n (t 1 ,t 2 , ...,t n ) = (z±(ti)z±(t 2 ) ■ ■ ■ z±(t n )) , h>t 2 >...>t n . 

It is convenient to use a recursive formula 

M n (t u t 2 , ...,t n ) = e- 2 "'^-^M n _ 2 (t 3 , ...,t n ). (1.33) 

which follows directly from Eqs. (1.31) and (1.30). Combining Eq. (1.33) with the Taylor expansion of 
Eq. (1.32) we obtain an exact integral equation for ip±[(5, t\: 



V>±(/M) = 1 ± i(v/2) / dtxe'^Pih) 
Jo 

-(v 2 /A) [ dh [ dhe-^-^PihWhy^frh}. 



(1.34) 



One can evaluate second time-derivative of both sides of Eq. (1.34) and transform this integral equation 
into a second order differential equation [25] 



dt 2 



+ 



2 7 



with initial conditions 



^±(0) = 1 



dln(3(t) 
dt 

di>± 



di/;± v 2 







dt 



t=o 



(1.35) 



(1.36) 



Dephasing of a qubit by 1 / f -noise 9 

2.2 Generating functional 

In the limit E{ <C T the random functions can be expressed through z + (t) or by z-(t) with 
equal probabilities. Using Eq. (1.17) we thus can rewrite the memory functional (1.27) in terms of the 
functionals ip±, Eq. (1.32): 

mt\ = lJ2 eTt(v/2) /ot m dt '^ ^(*). *] • (i-37) 
1 ± 

From Eqs. (1.37,1.35) follows the differential equation for i/j[f3, t]: 



d 2 ^ 

with initial conditions 



+ 



dlnP(t) . 
2 7 - - l vf3(t) 



— li_^ 7 ^ = o (1.38) 
at 



= */9(-0)£. (1.39) 

i=0 2 



Below we use this equation to analyze decay of the free induction and echo signals. Note that Eq. (1.38) 
is the Ei/T — > limit of equation derived in Ref. [25] for arbitrary Ei/T. 

For the free induction signal [3(t > 0) = 1 and /?(— 0) = 0. For this /^-function Eq. (1.38) with initial 
conditions (1.39) yields the following phase memory functional 

M) = ^ E E P f 1 - ^ + PI*) e-W*™* . (1-40) 

4 ^ g=±lp=±l V ^7 / 

where fi = 1 — (v/2'y) 2 . This solution is the i£j/T — > hmit of the result obtained in Ref. [17]. At short 
times, jt <C 1, Eq. (1.40) can be approximated as 

^ pm {t) = l-l[t . (1.41) 



However t/j prn given by Eq. (1.40) does not describe the free induction decay, exp iv f^(t')dt' 
even within our simplified model. We have to consider 

**(*>- »P (£)*t±^ (1.42) 

rather than tppm(t), because the latter neglects the dispersion of Eq, (1.17), which is due to the term 
u£(0). From Eqs. (1.35) and (1.36) it follows that 

^ fi (t) = e {iv/2 -^ ]t (^T 1 cosh 7/^ + sinh 7 ^) . (1.43) 

This expression is the Ei/T — > limit of the result obtained in Ref. [10]. Comparison of the free 
induction decay (1.43) with decay of the phase memory is presented in Fig. 1.2. This difference is 
especially important for the case of 1/ /-noise, when many fluctuators contribute. 

The calculation of the echo decay can be done in a similar way. The results for 2- and 3-pulse echo 
are, respectively (cf. with Ref. [17]): 

-27T12 

<M2r 12 ) = ^-nrE [(i + ^)(i±m)e ±2 ^ T12 

-(l-^)(lT ¥2 )e^] , (1.44) 

1 , , 2 



i>e?,{ri2 + Ti 3 ) = ^e2(2r 12 ) + - 7^— j- (COS 2/i 2 Ti 2 

-COS2/XXT12) fl - e- 27(Tl3 " Tl2)N ) e" 27T12 . (1.45) 



10 




Figure 1.2. Comparison of the phase memory function 
(1,1a) with the free induction signal (2,2a). v/2^y — 5 
(1,2), and 0.5 (2,2a). Time is measured in units of 

(ST" 1 )- 



Figure 1.3. Two-pulse echo signal for v/2~/ — 5 (1) 
and 0.8 (2). Time is measured in units of (2^y^ 1 ). 



Here T13 is the time between first and third pulse, \i\ + = yjl — (v/2j) 2 . The function ip e 2(2T\2) 
is shown in Fig. 1.3. Note that at v > 27 the time dependence of the echo signal shows steps similar to 
what was experimentally observed for the charge echo [2]. The expressions for the echo signal have a 
simple form at v S> 7. In particular, the expression for the two-pulse echo acquires the form 

<M2r 12 ) = e" 27T12 ([1 + (2-y/v) smvT 12 ] . 

Consequently, the plateaus occur at VT12 = kir/2 — arctan(27/i>) and their heights exponentially decay 
with the number k. 

3. Summation over many fluctuators 

To average over a set of the fluctuators we assume that dynamics of different fluctuators are not 
correlated, i.e., (£') = 0, unless i = j. Under this condition the generating functional is aproduct 

of the partial functionals, Hence, the generating functional can be presented as 

y[(3(t),t] = Y[^ {i) (t) = e £i = e -JC(t) _ (L46) 

i 

Since logarithm of the product is a self-averaged quantity, at large number of fluctuators, M 3> 1, one 
can replace the sum £\ In by J\f (In iP)f, where (■ ■ -)f denotes average over the fluctuators interlevel 
spacings Ei, their interaction strength, Vi, and tunneling parameters, sinr^j. Furthermore, for M >> 1 one 
can employ the Holtsmark procedure [31], i. e., to replace (IiiV')f by (ifi — 1)f, assuming that each of 
is close to 1. Thus, for large Af the approximate expression for KL(t) is 

K{t) ^(t)) = M (1 - m) F ■ (1-4V) 

i 

To evaluate K.(t) one has to specify the distributions of the parameters E^, v^, and Q; L that characterize 
the fluctuators. Taking into account only the fluctuators with E{ < T we can write the number of 
fluctuators per unit volume as PqT. It is natural to assume that the density of states, P is a T-independent 
constant as it is in structural glasses [14]. 

The conventional estimation of the distribution function of relaxation rates is based on the fact that 
the tunneling splitting A depends exponentially on the distance in real space between the positions of the 
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two-state fluctuate* (on the distance between the charge trap and the gate), as well as on the height of the 
barrier between the two states Assuming that parameters like distances and barrier heights are distributed 
uniformly, one concludes that A-distribution is oc A -1 . Since 7 is parameterized according to Eq. (1.20), 
this implies V{9) = 1/ sin 9. As a result, cf. with Ref. [14], 

V(E,9) = P /sm0. (1.48) 

The coupling constants, Vi, determined by Eq. (1.13), contain cos 0£ and thus are statistically correlated 
with 0j. It is convenient to introduce an uncorrected random coupling parameter, ui as 

Ui = g(ri)A(ni) , Vi =Ui cos 0* . (1.49) 

It is safe to assume that direction, rij, of a fluctuator is correlated neither with its distance from the 
qubit, rj, nor with the tunneling parameter 0j and thus to replace A(n) by its angle average, (|A(n)|) n . It 
is likely that the coupling decays as power of the distance, r: g(r) oc g/r b . If the fluctuators are located 
near a d-dimensional surface, then 

TV a\ -d/b-1 9 a d en\ 

™ = M U >* = ^>* = Tjmv*- (L50) 

Here ad is a dimensionless constant depending on the dimensionality d while is a typical distance 
between the fluctuators with E\ < T . In the following we will for simplicity assume that 

' min < r T < ' max i 

(1.51) 

where r m - m (r max ) are distances between the qubit and the nearest (most remote) fluctuator. Under this 
condition rj oc T b / d is the typical constant of the qubit - fluctuator coupling. As soon as the inequality 
(1.5 1) is violated the decoherence starts to depend explicitly on either r m j n or r max , i. e., become sensitive 
to mesoscopic details of the device. This case will be analyzed elsewhere. In the following we assume 
that d = b, as it is for charged traps located near the gate electrode, see Fig. 1.1. 

The dependences of the generating functional ip(/3, t\u, 7) on the coupling constants u and transition 
rates 7 of the fluctuators are determined by Eq. (1.27). Substituting equations (1.20) and (1.49) into 
(1.27) and the result - into (1.47) one obtains K(t) in the form: 

JC{t) = T 1 [ % I"' 2 —A l-^[/3,t|u cos 0,70 sin 2 9] } . (1.52) 
J u z J sm 9 L J 

This expression allows one to calculate the dephasing rate in the case of many surrounding fluctuators 
for various qubit manipulations. 

To start with consider the phase memory functional with (3(t') = 0(i'), i. e., free induction signal will 
£i(0) = 0. At small times, t 7Q" 1 , one can use Eq. (1.41) for ippm(t). The integration over u yields 



f°° du , , . . , 1 , q 
/ -j 1-Vpm (*) = 7 t 2 7osm 2 0. 
Jo u 2 4 



Performing the following integral over 9 one obtains /C pm oc rj^nt 2 . 

It is slightly trickier to estimate JC pm (t) at large times, 7ot 3> 1. One can show that the decoherence in 
this limit is due to the fluctuators, which coupling with the qubit is atypically weak: u ~ t~ x <C 7o- As a 
result, in the leading in 1/ (70*) approximation i/j prn (t) = cos ut/2. This asymptotics can be interpreted 
in the following way. At t 3> I/70 a typical fluctuator had flipped many times and its contribution to the 
qubit phase, which is proportional to 



J*m')-m]dt'\ oct : 
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and does not depend on 7 and, hence on 8. Therefore the integral over 9 in (1.52) diverges logarithmi- 
cally. The proper cut-off is determined by the condition jt « 1, i. e., is the value of 9 = 9 m i n {t), which 
allows approximately one flip during the time t. Using Eq. (1.20) we estimate 9 m i n (t) as (7ot) -1 / 2 ). 
This cut-off reflects the fact that the fluctuators with too low tunneling rates do not change their states 
during the measurement time t. 

The estimate for K, pm {t) can be then summarized as (cf. with Ref. [17]), 

!C pm (t)*v{ if + f 7 °^! ; d-53) 
pmw ' I fm7 i for 7 t > 1 . 

Now we can define the dephasing time r v by the condition KL{t^) = 1. Using Eq. (1.53), we get 

t v = max {tT 1 ln~ 1 (7o/??), (mo)^ 2 } ■ (1-54) 

The echo decay can be calculated in a similar way, cf. with Ref. [17]: 

£ e2 (2Ti 2 ) ~ W ■ mm{l,7 Ti 2 } , (1.55) 
£ e 3(Ti2 + Ti3) ~ 7]t u min{l , In T13/T12} . (1.56) 

The dephasing time for the two-pulse echo decay is then 

^ = max{r ? - 1 ,(r ? 7o)~ 1/2 } . (1.57) 

Let us discuss the physical meaning of the results (1.53)— (1.57). If there is no flips of the fluctua- 
tors, then the contribution of a given fluctuator to the total phase gain during the observation time is 
£(0) Jq P(t') dt'. During the time interval t <C 7^ each fluctuator can flip only once. If it flips at time 
t\ , the accumulated relative phase is ±2 f3(t') dt'. To obtain the total phase gain one has to average 
over all possible moments of flips: 



\6<p(t)\ ~7o / dh 



0(t')dt' 

Since v oc r -3 , nearest neighbors are important, and typical value of v is 77. Thus we immediately obtain 
K{t) ~ f?7o£ 2 - It can be shown [13] that in this case the random process is Markovian. Consequently 
in this case the situation can be characterized by a pair conditional probability K(E,t\Eo) to find the 
spacing E at time t under condition that at t = it was Eq. It has the Lorentzian form, 

K(E,t\E )= 1 - w -^ Tfm , (1.58) 

where T(x) ~ rj^fot. This time dependence can be easily understood in the following way. The nearest 
region of r, where at least one fluctuator flips during the time interval t gives the maximal contribution. 
The size of this region can be estimated from the condition i-b2V 3 7o^ « 1 that yields r~ 3 « PoTjot. 
The corresponding change in the qubit's interlevel spacing is then given by the interaction strength at this 
distance, 7770^. 

During a time interval t 3> 7 ( ^ 1 a substantial contribution comes from the fluctuators with 7^ <C 
7 _1 <C t, which experience many flip-flops. Since a fluctuator having £ = 1 can flip only to the state 
with £ = the contributions of successive hops are not statistically independent. The density of most 
important fluctuators is of the order Po^V 3m 7o^ and the substantial region of r is determined by the 
relation Po^ 3m 7o* ~ 1- As a result T(t) ~ r?ln7ot. This dependence holds until t < 7"^, where 
7min is the minimal relaxation rate in the system. At t > 7 mi 1 n the quantity T(t) saturates at the value 
of the order of Too ~ r\ ln(7o/7 m i n )- The random process in this case is non-Markovian and cannot be 
fully characterized by the pair correlation function (1.58). 
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The above estimates do not describe decay of the free induction signal due to beats between contribu- 
tions of different fluctuators. For t <C 7 ^ 1 this decay can be evaluated in the same way as influence of 
the static inhomogeneous broadening [12, 17]. Averaging Eq. (1.26) we get Kf(t) = T^t. At large time, 
t S> 7Q -1 , one can assume that the probability to find a fluctuator in a given state is just the equilibrium 
one. Thus 



Averaging this solution over the hopping rates and positions of the fluctuators we arrive at the same 
result. Consequently, the free induction signal decays much more rapidly than the phase memory. 

The important point is that at large observation time, t 3> 7q~ 1 , there are optimal fluctuators respon- 
sible for decoherence. The distance r op t(T), between the optimal fluctuators and the qubit is determined 
by the condition 



If the coupling decays as l/r b and the fluctuators are distributed in a d-dimensional space, then r d 1 dr — > 

V(v) oc V 

From this one concludes that at d < b the decoherence is controlled by optimal fluctu- 
ators located at the distance r opt provided they exist. At d > b the decoherence at large time is determined 
by most remote fluctuators with v = v m - m . If d < b, but the closest fluctuator has v max <C 70, then it is 
the quantity t> max that determines the decoherence. In both last cases JC(t) oc t 2 , and one can apply the 
results of Ref. [10], substituting for v either t> m i n or v max . Since r opt depends on the temperature there 
can exist a specific mesoscopic behavior of the decoherence rate. A similar mesoscopic behavior of the 
decoherence has been discussed for a microwave-irradiated Andreev interferometer [9]. 

4. Simulations 

The procedure outlined above leaves several questions unanswered. First, how many experimental 
runs one needs to obtain the ensemble-averaged result for a single fluctuator? Second, when the contri- 
butions of several fluctuators can be described by averages over the fluctuators' parameters? 

The second question is the most delicate. The situation with a qubit interacting with environment in 
fact differs from that of a resonant two-level system in spin or phonon echo experiments. In the first case 
the experiment is conducted using a single qubit surrounded by a set of fluctuators with fixed locations, 
while in the second case many resonant TLSs participate the absorption. Consequently, one can assume 
that each TLS has its own environment and calculate the properties averaged over positions and transition 
rates of the surrounding fluctuators. How many surrounding fluctuators one needs to replace the set of 
fluctuators with fixed locations (and transition rates) by an averaged fluctuating medium? 

4.1 One fluctuator 

To answer these questions we have made a series of simulations. First, using Eq. (1.27) we have cal- 
culated the two-pulse echo decay for a qubit with a single neighboring fluctuator. The random switching 
between the £ = and £ = 1 states of this fluctuator was modeled by the Poisson process with time 
constant 7. The result of simulations for u/27 = 5 are shown in Fig. 1.4. The left panel presents the 
average over 100 random runs, while the right one shows the average over 1000 runs. Ultimate averaging 
over infinite number of runs would give the analytical result (1.44) shown by the curve 1 in Fig. 1.3. One 
can see that averaging over 1000 runs is sufficient to reproduce the analytical result with good accuracy, 
in particular, to observe the characteristic plateau around. Note that the plateau is qualitatively similar to 
experimentally observed for the charge echo in Josephson qubits. [2] Note also that with only 100 runs 
made, the dispersion of the signal indicated by the error bars is huge. Therefore, experimentally one 
needs at least many hundred runs to obtain reliable averages. 




v(r opt ) « 70 (T) . 



(1-59) 
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Figure 1.4. Two-pulse echo decay (in logarithmic scale) for the case of a single neighboring fluctuator, u/27 = 5. The left 
panel corresponds to 100 realizations of the random switching process, while the right one corresponds to 1000 realizations. 

4.2 Check of the Holtsmark procedure 

To check the validity of the summation over different fluctuators using the Holtsmark procedure, we 
perform simulations for many fluctuators. The fluctuators are assumed to be uniformly distributed in 
space at distances smaller than some r max . Then, the normalized distribution function of the coupling 
constants and relaxation rates, V(u, 9) can be specified as 

V(u,6) = (Mu min /u 2 ) [sm6 ln(tan0 min /2)]' 1 (1.60) 

with u € {limin, co}, while 9 6 {O m in, 7r /2}- Here small 9 correspond to slow fluctuators, 7 = 70 sin 2 9. 
The quantity r) (1.50) that characterizes the fluctuator density is given now as 77 = Afiimin/ hi(tan 6> mm /2). 

The results of simulations for small times are shown in Fig. 1.5. For simplicity, in these simulations 
the transition rate was assumed to be the same, 7 = 70, for all fluctuators. In order to check the analytical 
result (1.53) for the free induction signal, it is convenient to plot K.(t)/rjt versus 2ft. One can see that 
the predicted asymptotic behavior works rather well for ft <C 1. By substituting (1.44) into (1.52) it 
can be easily shown that the two-pulse echo signal has a similar asymptotic for small t, K oc t 2 , only 
the coefficient is twice that for the free induction signal. This result is also perfectly reproduced by the 
simulations, which justifies the use of the Holtsmark procedure for ft -C 1, i. e. when K is small. The 
results for large times are shown in Fig. 1.6. Here it was important to take into account scatter in values 
of 7 because behavior at large t is controlled by numerous fluctuators that flip very seldom, i. e. have 
small 7. One observes that analytically predicted behavior of the phase memory for the free induction 
signal, K oc t In t at t S> 7(7 1 is fully confirmed by the simulations. The analytical result was obtained by 
making a rough cutoff of slow fluctuators at 9 = (to*) -1 / 2 that led to K, pm (t) « rjUnjot, see Eq. 1.53. 
From simulations we can see that a more accurate expression at large times is K, pm {t) m {qt/2) In 270?;, 
which corresponds to the straight line in Fig. 1.6. 



4.3 Many fluctuators with fixed locations 

The curves presented in Figs. 1.5 and 1.6 were calculated by averaging over many random sets of 
fluctuators. This allowed us to make a reliable check of the analytical results based on the Holtsmark 
procedure. The next step is to check whether it is appropriate to average over the fluctuators positions 
though in a real system the fluctuators' parameters are fixed. For this purpose we compare the results for 
three different sets of fixed fluctuators with different coupling constants distributed again according to 
Eq. (1.60), however with fixed 9. For each set of fluctuators, we find the average signal over 1000 runs 
and show it in Fig. 1.7. The results strongly depend on the fluctuator density 77 = Wu m i n - For 77 = 300, 
the right panel, different sets of fluctuators lead to similar behavior of K(t), rather close to the "expected" 
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-□ — two-pulse echo 
-• — free induction 





100 



Figure 1.5. Short-time phase memory decay of two- 
pulse echo (1) and free induction (2) for u m in = 27. 
Lines (3) and (4), correspond to analytically calculated 
slopes: 7r/4 and tt/8, respectively. The results are aver- 
aged over 25000 realizations of random telegraph noise 
in N — 10 fluctuators randomly distributed in space. 



Figure 1.6. Long-time phase memory decay of the free 
induction for u m in = O.O27. Solid line corresponds to 
(1/2) ln(27o£). The results are averaged over 5000 real- 
izations of random telegraph noise in N — 10 fluctuators 
randomly distributed in space. 
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Figure 1.7. Two-pulse echo decay for three different fluctuators sets Umm = 27, N = 10 (left panel), and 300 (right panel). 



behavior obtained by averaging over 5000 different sets. For larger densities the reproducibility is even 
better. However for rj = 10, the left panel, each set of fluctuators is characterized by a very specific 
type of the signal. The function K(t) for a given set usually differs dramatically from the "expected" 
behavior that we obtained by averaging over 5000 sets. For such a low fluctuator density, it is hopeless 
to fit the experimental data with our analytical formulas obtained by averaging over fluctuator ensemble, 
like Eq.(1.53). Note that small error bars on the plot mean only that the measured signal is reproducible if 
it is averaged over 1000 experimental runs. However, even a slightly different arrangement of fluctuators 
in the gate can produce a very different signal 

If the fluctuator density is low, or, in other words if we are in the mesoscopic limit, the signal should be 
essentially determined by one, the most important fluctuator. To check this we have calculated the echo 
signal in the presence of J\f fluctuators, and compared it with the signal in the absence of the strongest 
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Figure 1.8. Two-pulse echo signal for the case of a TV = 10 neighboring fluctuators with the same transition rate , 7. It is 
assumed that u m i n = 7. 

fluctuator. The results are shown in Fig. 1.8. One can see that one needs really many fluctuators to avoid 
strong mesoscopic fluctuations. More detailed studies of mesoscopic fluctuations are planned for future. 

5. Comparison with the noise in the random frequency deviation 

The conventional way is to express the environment-induced decoherence through the noise spectrum, 
S x {uj) = 2 $™ dte iwt ({X(t)x[o))() F . Using Eq. (1.17) we get 

Let us start averaging over fluctuators by integration over 0i. Since we are interested in small frequen- 
cies, we can replace sin# — ► 9, cos 6 — ► 1 and replace the upper limit tt/2 of the integration by infinity. 
In this way we get Sx{uj) = ir , qu TQ3X /2u:. Here we have taken into account that the summation over 
the fluctuator strength, u, is divergent at the upper limit corresponding to the minimal distance, r m - m , 
between the fluctuator and the qubit. Thus, the closest fluctuators are most important. We observe that 
our model leads to 1// noise in the random force acting upon the qubit. However, the noise is mainly 
determined by the nearest fluctuators, while decoherence (at long times) is dominated by the fluctua- 
tors at the distance r opt given by Eq. (1.59). Since r opt S> r m - m , the decoherence cannot, in general, be 
expressed only through Sx(oS). 

Now we can compare our result given by Eq. ( 1.53). From Eq. (1.5) one obtains (cf. with Ref. [5]) 

K{t) = (l/2)r/u max t 2 | \uu)i r t\ 

where Wj r is the so-called intrinsic infrared cutoff frequency for the l/f noise. [3] It is clear that the 
results differ significantly. Even in the case when 7oi <C 1 when the random process is Markovian, 
results differ both by order of magnitude and temperature dependence. The reason for this discrepancy 
is that dephasing and l//-noise are determined by different sets of fluctuators. 

6. Applicability range of the model 

Let us discuss the applicability range for the used approach. Firstly, random fields acting on the qubit 
were assumed to be classical. This is correct provided the typical hopping rate of a fluctuator, 70, is much 
less than its typical interlevel spacing, E ~ T. This is the case, indeed, for fluctuator interaction with 
both with phonons and electrons because of weak coupling. Secondly, we did not discuss the mechanism 
of decoherence due to direct translation of excitation from the qubit to fluctuators. The most important 
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part of such interaction can be written as 

i 

where a± = a x ± icr^, and r± = ± ir y . One can expect that the coupling constant Ui is of the same 
order of magnitude as Uj, i. e. ~ g/r 3 . Assuming the constant fluctuator density of states P$ we can 
estimate the typical energy defect for translation of an excitation to the fluctuator at the distance r as 
(6E) ~ (iir 3 ) -1 . The effect of the above off-diagonal interaction can be then estimated as U/(5E) « 
Pod ~ t]/T. This ratio should be small within the applicability range of our theory. Indeed, qubit will not 
be useful if its characteristic decay rate, 77, exceeds its interlevel spacing Eq, which should be, in 
turn, much less than the temperature. However, at very low temperatures and at not too long observation 
times the above processes could be important [32]. 

Another issue, which has not been analyzed, it the decoherence near the degeneracy points where 
cos 9 q — > 0. These points are of specific importance since linear coupling between the qubit and the 
fluctuators vanishes. The conventional way, see, e. g., Ref. [33]) is to introduce the model coupling as 
V2 = XX 2 (t) a z . We believe that the model still needs a careful derivation. 

7. Conclusions 

The simple model discussed above reproduces essential features of the decoherence of a qubit by 
"slow" dynamical defects - fluctuators. This model is valid for qubits of different types. The main 
physical picture is very similar to that of the spectral diffusion in glasses. 

The phase memory decay is due to irreversible processes in the fluctuator system. In the case on 
ensemble-averaged measurements it can be directly determined from the echo-type experiments. In the 
experiments of the free-induction type the decay of the signal is due both to the finite phase memory and 
to the beats between different values of the qubit eigenfrequency. 

The effective rate of the phase memory decay depends on the relation between the typical interaction 
strength, r/(T), and the typical fluctuator relaxation rate 70 (T). The first quantity is just a typical devi- 
ation of the qubit eigenfrequency produced by a thermal fluctuator (with the inter-level spacing of the 
order of temperature T) located at a typical distance. The second quantity is the maximal flip-flop rate 
for the thermal fluctuators. 

The estimates for the phase memory time in the limiting cases 77 3> 70 and 77 <C 70 are summarized in 
Eq. (1.54). In the first limiting case during the decoherence time only few fluctuators flip. Consequently, 
the decoherence is governed by Markovian processes. In the opposite limiting case, typical fluctuators 
experience many flip-flops during the decoherence time, the subsequent deviations in the qubit frequency 
being statistically -dependent. As a result, the decoherence process is essentially non-Markovian. 

The details of the decoherence depend strongly on the concrete type of the fluctuators, namely on the 
distribution of their flip-flop rates, on the range of their effective field acting upon the qubit, and on the 
distribution of the fluctuators in real space. 
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